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(Dated: February 7, 2008) 

Branching fraction and asymmetry measurements of charmless B + — > K'^h^h^ (where /ii,2 = 
K, 7r) decays are presented, using a data sample of 232 million T(4S) — > BB decays collected with 
the BABAR detector at the SLAC PEP-II asymmetric-energy B factory. Using a maximum likelihood 
fit, the following branching fraction results were obtained: B{B + — > K* + K + K ~) = (36.2 ± 3.3 ± 
3.6) x 10~ B and B(B + -> K'+tt+ty') = (75.3 ± 6.0 ± 8.1) x 10" 6 . Upper limits were set for B{B+ 
-> K* + ir + K~) < 11.8 x 10~ 6 and B{B + -> K'+K+ty') < 6.1 x 10~ 6 at 90% confidence level. 
The charge asymmetries for the decays B + — > K* + K + K~ and B + — » K* + it + -k~ were measured to 
be Ak'kk = 0.11±0.08±0.03 and Art*™ = 0.07± 0.07± 0.04, respectively. The first error quoted 
on branching fraction and asymmetry measurements is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charmless decays of B mesons to three-body final 
states are very important in aiding the understanding 
of the weak interaction and complex quark couplings de- 
scribed by the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix elements Q. Improved experimental measurements 
of these charmless decays, combined with theoretical de- 
velopments, can provide significant constraints on the 
CKM matrix parameters or uncover evidence for physics 
beyond the Standard Model. For example, the branch- 
ing fraction of the decay B + — > K* + ir + K~ is sensitive 
to CKM matrix elements V t d and V u b (see Figure^). Ad- 
ditionally, a B + — ► K* + n + K~~ branching fraction value 
equal to or smaller than the branching fraction of the 
Standard Model suppressed decay B + — ► K* + K + w~ 
would be an indication of new physics. 




FIG. 1: Penguin (a) and tree (b) Feynman diagrams for the 
decay B + -» K* + ir + K~. 

We identify K* + mesons through their decay to K®ir + . 
Charged B-meson decays to K^tt + (where hi 2 = K 

or tt) arc dominated by -ftT* + /i^~h^~, but can also proceed 
via a nonresonant component as well as through observed 
intermediate charmless resonances such as B + — > K* + <p 
or B + — > K* + p° [2J, Li |j|, or other as- yet-unobserved 
intermediate charmless resonances. To date, there are 
only limits on the charmless decays B + — > K* + h^h2 , 
measured by the ARGUS experiment 5] using less than 
0.2 fb -1 . 

Asymmetry measurements of charmless B decays can 
be used to probe for CP violation where the CP asym- 



metry is defined as: 



1 K*-h~h+ ' 1 K-> + h+h~ 

and r is the partial width of the charged B decay in the 
subscript. 

In the analyses presented in this paper, branching frac- 
tions of B+ -> K*+K+K- and B+ -> K*+tt+t:~ were 
measured for the first time, and upper limits were set for 
B + — ► K* + n + K~ and the Standard Model suppressed 
decay B + — > K* + K + ix~ , where charge-conjugate decays 
are also implied. The selection criteria required events 
with a reconstructed K^^h^h^ final state such that 
the total charmless contribution to the K* + h^Ii2 Dalitz 
plot could be measured (with charmed and charmonium 
resonances removed), including contributions from reso- 
nant charmless substructure. Finally, the AK*hth 2 va l - 
ues for the observed decays B + — * K* + K + K~ and 
B + —t K* + it + tt~ were measured. 

The data used in this analysis were collected at 
the PEP-II asymmetric-energy e + e~ storage ring with 
the BABAR detector 6]. The BABAR detector consists 
of a double-sided five-layer silicon tracker, a 40-layer 
drift chamber, a Cherenkov detector, an electromagnetic 
calorimeter and a magnet with instrumented flux return. 
The data sample has an integrated luminosity of 210 fb -1 
collected at the T(4S) resonance, which corresponds to 
(231.8 ± 2.5) x 10 6 BB pairs. It was assumed that the 
T(4S) decayed equally to neutral and charged B-meson 
pairs. In addition, 21.6 fb -1 of data collected at 40 MeV 
below the T(4S) resonance were used for background 
studies. 

Candidate B mesons were reconstructed from three 
tracks and a K®, where the K® was reconstructed from 
7r + 7r~ candidates. Each of the three tracks that were not 
associated with the K® were required to have at least 
12 hits in the drift chamber, a transverse momentum 
greater than 100 MeV/c and to be consistent with orig- 
inating from the beam-spot. These tracks were identi- 
fied as either pion or kaon candidates using energy loss 
(dE/dx) measured in the tracking system and the num- 
ber of photons measured by the Cherenkov detector and 
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their corresponding Cherenkov angles. Furthermore, the 
tracks were required to fail the electron selection based on 
dE /Ax information, their ratio of energy in the calorime- 
ter to momentum in the drift chamber, and the shower 
shape of the signal in the calorimeter. The K® candi- 
dates were required to have a reconstructed mass within 
15 MeV/c 2 of the nominal K° mass JtJ, a decay ver- 
tex separated from the B + decay vertex by at least five 
standard deviations, and a cosine of the angle between 
the line joining the B and K® decay vertices and the K® 
momentum greater than 0.999. 

To characterize signal events, three kinematic vari- 
ables and one event-shape variable were used. The 
first kinematic variable AE, is the difference between 
the center-of-mass (CM) energy of the B-candidate 
and y/s/2, where y 7 ! is the total CM energy. The 
second is the beam-energy-substituted mass tties = 
yj (s/2 + p i • Pb) 2 /Ef — p 2 B , where ps is the B momen- 
tum and (2?j,pj) is the four-momentum of the T(4S) in 
the laboratory frame. The third kinematic variable is 
the K°ir + invariant mass, rriK" > used to identify K* + 
candidates. Using these three kinematic variables, candi- 
dates were required to be in the ranges \AE\ < 0.1 GeV, 
5.25 < m ES < 5.29GeV/c 2 and 0.772 < m K . < 0.992 
GeV/c 2 . The event-shape variable is a Fisher discrimi- 
nant (J-) , constructed from a linear combination of the 
cosine of the angle between the B-candidate momentum 
and the beam axis, the cosine of the angle between the 
B-candidate daughters thrust axis and the beam axis, 
and the zeroth and second angular moments of energy 
flow about the thrust axis of the reconstructed B. 

Continuum quark production (e + e~ — > qq where q = 
u,d,s,c) was the dominant source of background. This 
was suppressed using another event-shape variable which 
was the cosine of the angle 9t between the thrust axis 
of the selected B-candidate and the thrust axis of the 
rest of the event. For continuum background, the distri- 
bution | cos 9t I is strongly peaked towards unity whereas 
the distribution is flat for signal events. Therefore, the 
relative amount of continuum background was reduced 
by requiring | cos#t| < 0.8. 

Simulated Monte Carlo (MC) events were used to 
study background from other B-meson decays. The 
largest B-background for B + — > K* + tt + it~ candidates 
comes from decays including charmonium mesons such 
as J/ipK*+, XcoK* + and ip(2S)K*+, where the char- 
monium meson decays to which are misidenti- 
fied as pions, or where the charmonium meson decays 
directly to 7r + 7r~. These background events were re- 
moved by vetoing reconstructed 7r + 7r~ masses in the 
range 3.04 < m„+ ir - < 3.17 GeV/c 2 , 3.32 < m ie + ie - < 
3.53 GeV/c 2 and 3.60 < m n +„- < 3.78 GeV/c 2 , iden- 
tifying the J/ip, Xco and ip(2S) mesons, respectively. 
For B+ -> K*+K+K- candidates, J/^K+ and XcoK*+ 
events were removed by rejecting events with a recon- 
structed invariant mass in the range 3.04 < m K *+ K - < 



3.17 GeV/c 2 and 3.32 < m K+K - < 3.53 GeV/c 2 , respec- 
tively. 

Potential charm contributions from B + — > D Q (— > 
K* + h^)h2 events were removed from corresponding 
B + — > K* + fi2 hi candidates by vetoing events with 
a reconstructed K* + hi invariant mass in the range 
1.83 < niK'h < 1.91 GeV/c 2 . Additional decays such 
as B+ -> D~(-> #°7r-)7r+7r+ and B+ -> AT*+B°(-> 
K + ir~) were removed from B + — * K* + tt + tt~ and B + — > 
K* + tt + K~ I K* + K + ir~ candidates, respectively by veto- 
ing events with reconstructed Kit invariant masses in the 
range 1.83 < < 1.91 GeV/c 2 , using the same veto 

range for D° and D + candidates. Studies of MC events 
showed that the largest remaining charmed background 
was B+ -> D°(-> K* + tt-)it + , with 18% of these events 
passing the veto. Surviving charmed events had a recon- 
structed D mass outside the veto as a result of using the 
wrong 7T + , K + or which was incorrectly selected from 
the other B decay in the event. 

After the above selection criteria were applied, a frac- 
tion of events for all decays had more than one candi- 
date. For those events, one candidate alone was selected 
by choosing the candidate whose fitted B decay vertex 
had the smallest \ 2 value. Studies of MC events showed 
the selection of B + — > K* + tt + tt~ events produced the 
largest number of multiple candidates, in 29% of events, 
where for these multiple candidates the correct one was 
reconstructed 70% of the time. 

After all requirements, there were five main sources of 
-B-background: two-body decays proceeding via a char- 
monium meson; two and three-body decays proceeding 
via a D meson; combinatorial background from three 
unrelated particles (i^* + /i^"/i^"); charmless two or four- 
body B decays with an extra or missing particle and 
three-body decays with one or more particles misidenti- 
fied. Along with selection efficiencies obtained from MC 
simulation, existing branching fractions for these modes 
|9j were used to estimate their background contributions 
which were included in fits to data. 

In order to extract the signal event yield for the channel 
under study, an unbinned extended maximum likelihood 
fit was used. The likelihood function for N candidates is: 

^ / m \ n / m \ 
c = ah exp y ^ ni J n \ Yl n ^ P ^ d ^ £ o)j . ( 2 ) 

where i and j are integers, M is the number of hypothe- 
ses (signal, continuum background and B-background), 
rii is the number of events for each hypothesis determined 
by maximizing the likelihood function and Pi(a,Xj) is 
a probability density function (PDF) with the param- 
eters a associated with x, where x can be any of the 
four variables tues, AE, J 7 , and tuk*- The PDF is 
a product P{(a,&) = Pi(a mES ,m ES ) ■ Pi(aAE,AE) ■ 
Pi{oijr,!F) ■ Pi(a mK , ,rnif»). Studies of MC simulations 
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showed correlations between these variables were small 
for signal and continuum background hypotheses. How- 
ever for B-background, correlations were observed be- 
tween toes an d AE, which were taken into account by 
forming a 2-dimensional PDF for these variables. The 
parameters of the signal and i3-background PDFs were 
determined from MC simulation. The continuum back- 
ground parameters were allowed to vary in the fit, to 
help reduce systematic effects from this dominant event 
type. Upper sideband data, defined to be in the region 
0.1 < AE < 0.3 GeV and 5.25 < to es < 5.29GeV/c 2 , 
was used to model the continuum background PDFs. 
For the toes PDFs, a Gaussian distribution was used for 
signal, a threshold function 10] for continuum and the 
combination of a Gaussian and threshold function for B- 
background. For the AE PDFs, a sum of two Gaussian 
distributions with the same mean was used for the signal, 
a first-order polynomial for the continuum background 
and the combination of the two Gaussians and a first- 
order polynomial was used for £?-background. The T sig- 
nal, continuum and _B-background PDFs were described 
using the sum of two Gaussian distributions with distinct 
means and widths. Finally for rriK* PDFs, the sum of 
a Breit-Wigner and a first-order polynomial was used to 
describe the signal, continuum and B-background distri- 
butions. The first-order polynomial component of the 
rriK' PDFs was used to model misreconstructed events 
for signal and background. 

In order to allow uncertainties and corrections due to 
MC simulation to be calculated and applied to the signal 
modes under study, the decay B + — ► > K* + n~)ir + 
was used as a calibration channel. These events were 
selected using the B + — > K* + ir + n~ selection criteria, 
but requiring the reconstructed K* + tt~ invariant mass 
be in the range 1.84 < m K , +1T - < 1.88GeV/c 2 . With 
more than 1800 signal events and approximately a 4 to 1 
signal to background ratio in the total number of B + — > 
D°{— * K* + h~)tt + candidates, it was possible to fit the 
signal PDF parameters for this mode. 

Branching fractions, B, are usually calculated using 
the following equation, B = n S i g /{N BB - x e), where n S i g 
is the fitted number of signal events, e is the average 
signal efficiency obtained from MC simulation and N BB 
is the total number of BB events. For the charmless 
B + — > K* + h^h2 branching fraction, the average effi- 
ciency cannot be taken directly from MC events. This 
was because the efficiency varies over the Dalitz plane 
and the distribution of events in the Dalitz plane is un- 
known before fits to data. To calculate the branching 
fraction, a weight was assigned to each event, j, as: 

where k is an integer and V S ig,i is the signal row of the 
covariance matrix obtained from the fit |ll| . This proce- 



dure is effectively a background subtraction where these 
weights have the property • VVj = n S i g . The branching 
fraction is then calculated as B — Wj/(ej x N BB ) 
where e., (a function of mr r . , ,_, to 2 , + and the K* + — > 

J v K* + h 1 ' h 1 hj 

K®ir + decay helicity angle, 7~Ik*+) varies across phase 
space and is simulated in bins using over seven million 
MC events for each channel. The size of the bins are op- 
timized to provide continuous coverage of the efficiency 
distribution. 

Figure [21 shows the fitted projections for B + — > 
K*+K+K-, B+ -> K*+ir+K-, B+ -> K*+K+n- and 
B + — > K* + ir + ir~ candidates, while the fitted signal 
yield, measured branching fractions, upper limits and 
asymmetries are shown in Table [I] The candidates in 
Figure |21 are signal-enhanced, with a requirement on the 
probability ratio V S i g /{V S i g + Vbkg), optimized to en- 
hance the visibility of potential signal, where V S i g and 
Vbkg are the signal and the total background probabil- 
ities, respectively (computed without using the variable 
plotted). The 90% confidence level (C.L.) branching frac- 
tion upper limits (Stjl) were determined by integrating 
the likelihood distribution (with systematics uncertain- 
ties included) as a function of the branching fraction from 
to B VL , so that J BuL CdB = 0.9 J °° CdB. 




FIG. 2: Maximum likelihood fit projections of rriES for signal- 
enhanced samples of charmless B + — » K'^h^li^ candidates. 
The dashed line is the fitted background PDF while the solid 
line is the sum of the signal and background PDFs. The 
points indicate the data. The plot labeled (a) shows a projec- 
tion of B + -> K*+K+K- candidates, (b) B+ -> K* + -k + K~ 
candidates, (c) B + — » K* + K tv~ candidates and (d) B + — > 
K* + n + ir~ candidates. 

Contributions to the branching fraction systematic er- 
ror are shown in Table [H] Errors due to charged track- 
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TABLE I: Signal yields, efficiencies and branching fractions for B + K* + K + K~, B + -» A"* + tt + A"~, B+ -» i^*+if + 7r" and 
B + —* K* + tt + it~ , measured using Kgir + hih^ events. The first error is statistical and in the case of the measured branching 
fractions the second error is systematic. These efficiencies have taken into account that B(K* + — > K°n + ) = 2/3, assuming 
isospin symmetry, as well as B(K° — > Kg) and B(Kg — > 7r + 7r~) Q. The branching fraction upper limits at a 90% C.L. are 
shown for B + — > i("* + 7r + A - and _B + — > K* + K + ir~ . Asymmetries are reported (first error is statistical and the second is 
systematic) only for the channels with statistically significant yields. 





Mode 




Signal Events Efficiency Measured Branching Fraction Upper Limit Asymmetry [A.K*hh) 
Yield (%) " (x 10" 6 ) (x 10" 6 ) 


B + 


- K* + K+ 


K~ 


288 ± 26 


3.4 


36.2 ± 3.3 ± 3.6 


0.11 ± 0.08 ± 0.03 


B+ 


-> K* + TV + 


K~ 


20.1 ± 24.7 


3.5 


2.5 ± 3.1 ± 5.3 


11.8 


B+ 


-> K* + K~ 


"n- 


9.7 ± 17.1 


3.5 


1.2 ± 2.1 ± 2.0 


6.1 


B+ 






583 ± 46 


3.3 


75.3 ± 6.0 ± 8.1 


0.07 ± 0.07 ± 0.04 



ing efficiency and K® reconstruction efficiency were as- 
signed by comparing control channels in MC events and 
data. The error in the efficiency was due to limited 
MC statistics, generated for each of the decays B + — > 
K*+K+K-, B+ -► K*+tt+K-, B+ -► K* + K + %~ and 
B + — > K* + tt + it~ . Using a sample of e + e~ — > /i + /J~ 
decays, the uncertainty in the number of BB events 
was calculated to be 1.1%. To calculate errors due to 
the fit procedure, a large number of MC samples were 
used, containing the amounts of signal and continuum 
events measured in data, and the estimated number of 
-B-background events. The differences between the gen- 
erated and fitted values using these samples were used 
to ascertain the sizes of any biases. Biases of +4.2, 
+10.7, +5.1 and +6.4% were observed in the fitted sig- 
nal yields of B+ -> K* + K + R-, B+ -> K* + ir + K~-, 
B+ -> K*+K+tt- and B+ -> K* + n + ir~, respectively, 
that were a consequence of small correlations between fit 
variables. These biases were applied as corrections to ob- 
tain the final signal yields and half of the correction was 
added as a systematic uncertainty. The uncertainty of 
the ^-background contribution to the fit was estimated 
by varying the known branching fractions within their 
errors. Each background was varied individually and the 
effect on the fitted signal yield was added in quadrature 
as a contribution to the uncertainty. For the my. fit 
range there was also the possibility of -B-background con- 
tributions from nonresonant and higher K* + resonances 
which was modeled in the data fit using a LASS param- 
eterization |l2t lb'{l| . The contribution from this back- 
ground was estimated by extrapolating a Kit invariant 
mass projection fitted in a higher-mass region (0.992 < 
rriK* < 1.6 GeV/c 2 ), into the signal region. This esti- 
mated background was modeled in the final data fit, and 
half of the contribution to the signal was added as a sys- 
tematic error contribution. The extrapolation assumed 
there were no integrated interference effects between the 
Ktt background and the if*" 1 " (892) signal. The uncer- 
tainty due to reconstructing the wrong B + — > K* + h\h 2 
signal candidate as a consequence of K/ir misidentifica- 



tion, for example B + — * K* + K + K~ events being re- 
constructed as B + — > K* + n + K~ candidates was deter- 
mined using MC events and added as a systematic. The 
uncertainty due to PDF modeling was estimated from 
the calibration channel B + — > D°(— > K* + ir~)ir + and 
by varying the PDFs according to the precision of the 
parameters obtained from the calibration channel fit to 
data. In order to take correlations between parameters 
into account, the full correlation matrix was used when 
varying parameters. All PDF parameters that were orig- 
inally fixed in the fit were then varied in turn, and each 
difference from the nominal fit was combined in quadra- 
ture and taken as a systematic contribution. 



TABLE II: Summary of systematic uncertainty contributions 
to the branching fraction measurements B + — > K'^h^h^ ■ 
Multiplicative errors are shown as a percentage of the branch- 
ing fraction and additive errors are shown in events. The fi- 
nal column shows the total systematic error on the branching 
fraction. 



Error K 


"+K+K- 


K*+tv+K- 


K*+K+tt- 




source error(%) 


error(%) 


error(%) 


error (%) 


Multiplicative errors (%) 








Tracking 


2.4 


2.4 


2.4 


2.4 


K% Efficiency 


1.2 


1.2 


1.2 


1.2 


Efficiency 


5.3 


9.4 


8.2 


5.6 


No. of BB 


1.1 


1.1 


1.1 


1.1 


Tot. mult.(%) 


6.0 


9.8 


8.7 


6.3 


Additive errors 


(events) 








Fit Bias 


6.0 


1.1 


0.3 


18.7 


B-background 


6.9 


33.4 


1.3 


37.9 


A(892)*+ bkg 


19.5 


12.0 


4.5 


25.6 


Signal mis-id 


0.0 


14.2 


15.5 


0.0 


PDF params. 


8.4 


19.4 


3.0 


14.0 


Tot. add. 










(events) 


23.1 


42.9 


16.5 


51.4 


Total (10~ b ) 


3.6 


5.3 


2.0 


8.1 


For the decay 


s B+ -> 


K*+hfh 2 


", interference effects 


between the K*~ 


f (892) and S-wave 


final states 


(nonreso- 
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nant and -ftTQ + (1430)) integrate to zero if the acceptance 
of the detector and analysis is uniform; the same is true 
of the interference between the K* + (892) and D-wave 
final states (JC] + (1430)). Studies of MC events showed 
the efficiency variations were small enough to make these 
interference effects insignificant. The integrated inter- 
ference between JT* + (892) and other P-wave amplitudes 
such as i<f* + (1410) is in principle nonzero, but in prac- 
tice is negligible due to the small branching fraction of 
if*+(1410) -> K° s tt+ (6.6 ± 1.3% 7J) and the fact that 
the Kit mass lineshapes have little overlap. 

The CP-violating charge asymmetries for the decays 
B + — > K* + K + K~ and B + — > K* + n + ir~ were measured 
to be Ak*kk = 0.11 ± 0.08 ± 0.03 and A K >™ = 0.07 ± 
0.07±0.04, respectively, where the first errors are statisti- 
cal and the second errors are systematic. The background 
asymmetries A^l g KK and A^* 9 ^! which were expected 
to be consistent with zero, were measured to be 0.00 ± 
0.02 and 0.01 ± 0.01, respectively. As a further study, 
the asymmetry Adtt for B + — > D°(— * iC* + 7r~)7r + , also 
expected to be consistent with zero, was measured to be 
—0.01 ± 0.02 with a corresponding background asymme- 
try A 1 ^ 9 of —0.01 ± 0.06 (statistical errors only). 

The systematic error on AK*hih 2 was calculated by 
considering contributions due to track finding, fit bi- 
ases, B-background uncertainties and particle interac- 
tion asymmetries. The error due to fit biases was found 
to be negligible. Uncertainties due to charged tracking 
efficiency were assigned by comparing control channels 
in MC simulation and data. The contribution from B- 
background was calculated by varying the number of 
expected events within errors and by assuming a CP- 
violating asymmetry of ± 0.2, as there are no available 
measurements for these decays. The CP asymmetry as- 
sumed for the individual £>-backgrounds was chosen to 
be greater than any asymmetry observed for charged B 
decays. The interaction asymmetry of matter and anti- 
matter with the detector was studied using background 
measurements A^» 9 KK and AgJ^ v , where no effect was 
observed and the statistical precision of the measurement 
was added as a systematic contribution. 

In summary, we have analyzed K®tt + final states 

to obtain the first branching fraction measurements for 



B+ -> K*+K+n-. 
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the decays £r 



K* + K+K~ and B+ 



K 



*+_+- 



and no evidence for new physics was found, placing upper 
limits on the branching fractions B + — » K * + n + K~ and 
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